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In the vast literature on tsunami research, few articles have been devoted to en- 
ergy issues. A theoretical investigation on the energy of waves generated by bottom 
motion is performed here. We start with the full incompressible Euler equations in 
the presence of a free surface and derive both dispersive and non-dispersive shallow- 
water equations with an energy equation. It is shown that dispersive effects only 
appear at higher order in the energy budget. Then we solve the Cauchy-Poisson 
problem of tsunami generation for the linearized water wave equations. Exchanges 
between potential and kinetic energies are clearly revealed. 
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1. Introduction 

Oceanic waves can be devastating as shown by recent events. Whilst some areas 
are more vulnerable than others, the recent history shows that catastrophic waves 
can hit even where they are not expected. The tsunami waves generated by the 
huge undersea earthquake in Indonesia on 26 December 2004 caused devastation 
across most of the coasts of the Bay of Bengal. The tsunami waves generated by the 
massive submarine landslide in Papua-New Guinea on 17 July 1998 as well as the 17 
July 2006 Java tsunami and the 2 April 2007 Solomon Islands tsunami also caused 
devastation, but on a smaller scale. Unfortunately, such cataclysmic tsunamis are 
likely to be generated again by earthquakes, massive landslides or volcano eruptions 
[Synolakis and Bernard, 2006]. 

Information on tsunami energy can be obtained by applying the normal mode 
representation of tsunami waves, as introduced by Ward [1980]. For example, Okal 
[2003] considers the total energy released into tsunami waves. He obtains expres- 
sions for the energy of tsunamis (see his expressions (31) for a tsunami generated 
by an earthquake and (36) for a tsunami generated by a landslide). In the case 
of a landslide, he computes the ratio between tsunami energy and total change in 
energy due to the slide. In the present paper, we use the incompressible fluid dy- 
namics equations. Tsunamis have traditionally been considered as non-dispersive 
long waves. However various types of data (bottom pressure records [Ritsema et al., 
1995]; satellite data [Kulikov et al., 2005]; hydrophone records [Okal et al., 2007]) 
indicate that tsunamis are made up of a very long dispersive wave train, especially 
when they have enough time to propagate. These waves travel across the ocean 
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surface in all directions away from the generation region. Recent numerical com- 
putations using dispersive wave models such as the Boussinesq equations show as 
much as 20% reduction of tsunami amplitude in certain locations due to dispersion 
[Dao and Tkalich, 2007] . But one has to be careful with the interpretation of satel- 
lite data: as indicated by Kanoglu and Synolakis [2006], the mid-ocean steepness 
of the 2004 Sumatra tsunami measured from satellite altimeter data was less than 
10~ 5 . Nonlinear dispersive theory is necessary only when examining steep gravity 
waves, which is not the case in deep water. 

The wavelength of tsunamis and, consequently, their period depend essentially 
on the source mechanism. If the tsunami is generated by a large and shallow earth- 
quake, its initial wavelength and period will be greater. On the other hand, if the 
tsunami is caused by a landslide (which happens less commonly but can be devas- 
tating as well), both its initial wavelength and period will be shorter as indicated 
for example by Kulikov et al. [1996]. From these empirical considerations one can 
conclude that dispersive effects are a priori more important for submarine landslide 
and slump scenarios than for tsunamigenic earthquakes. 

Once a tsunami has been generated, its energy is distributed throughout the wa- 
ter column. Clearly, the more water is displaced, the more energetic is the tsunami 
(compare for example the December 2004 and March 2005 Sumatra tsunamis). 
Due to the large scale of this natural phenomenon and limited power of comput- 
ers, tsunami wave modellers have to adopt simplified models which reduce a fully 
three-dimensional (3D) problem to a two-dimensional (2D) one. This approach is 
natural, since in the case of very long waves the water column moves as a whole. 
Consequently the flow is almost 2D. Among these models one can mention the non- 
linear shallow water equations (SWE), Boussinesq type models, the Green-Naghdi 
and Serre equations. There is a wide variety of models, depending on whether or 
not the effects of run-up/run-down, bottom friction, turbulence, Coriolis effects, 
tidal effects, etc, are included. 

Today scientists can easily predict when a tsunami will arrive at various places 
by knowing source characteristics and bathymetry data along the paths to those 
places. Unfortunately one does not know as much about the energy propagation 
of such waves. Obviously tsunami amplitude is enhanced over the major oceanic 
ridges. Titov et al. [2005] clearly describe the waveguide type effect from mid- 
ocean ridges that has funnelled the 2004 megatsunami away from the tip of Africa. 
As emphasized by Kowalik et al. [2007], travel-time computation based on the 
first arrival time may lead to errors in the prediction of tsunami arrival time as 
higher energy waves propagate slower along ridges. At the beginning, the energy 
is essentially potential, although it depends on the generation mechanism. Then it 
redistributes itself into half kinetic and half potential energies. Finally, it converts its 
potential component into kinetic energy. How do these conversions take place? The 
purpose of this study is to shed some light on this topic and to see if the importance 
of dispersion in tsunamis can be studied by looking at the energy rather than at 
wave profiles. 

Previous researchers have considered the topic of tsunami wave energy. Dotscnko 
and Korobkova [1997], Kajiura [1970], Velichko et al. [2002] studied the energy 
exchange between the solid bottom and the overlying water associated with the 
bottom deformation. There were recent attempts to obtain equations for tsunami 
energy propagation. We can mention here the work of Tinti and Bortolucci [2000] 
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devoted to idealized theoretical cases and the work of Kowalik et al. [2007] using 
the energy flux point of view to study the changes in the 2004 tsunami signal as it 
travelled from Indonesia to the Pacific Ocean. We believe that these models can be 
improved, given the present state of the art in wave modelling. 

A point of interest is that some of the equations used for wave modelling have an 
infinite number of conserved quantities. There has been some confusion in the liter- 
ature on which quantities can be called energy. Indeed there is here an interesting 
question. In incompressible fluid mechanics, the internal energy equation is decou- 
pled from the equation of continuity and from the fundamental law of dynamics. 
It is used only when one is interested in computing the temperature field once the 
velocity distribution is known. In addition to the internal energy equation, one can 
write a total energy (internal energy + kinetic energy) equation, or a total enthalpy 
equation. The confusing part is that for perfect fluids one usually defines the total 
energy differently: it is the sum of internal energy, kinetic energy, and potential 
energies associated to body forces such as gravitational forces and to the pressure 
field. If in addition the fluid is incompressible, then the internal energy remains 
constant. In the classical textbooks on water waves [Johnson, 1997, Stoker, 1958], 
one usually introduces the energy E as the sum of kinetic and potential energies 
and then looks for a partial differential equation giving the time derivative dE/dt 
(incidentally the meaning of d/dt is not always clearly defined). In any case, when 
one uses a depth- integrated model such as the nonlinear SWE, one can compute 
the energy a posteriori (the potential energy is based on the free-surface elevation 
and the kinetic energy on the horizontal velocity). But one can also apply the non- 
linear shallow water assumptions to the full energy equation to begin with. Then 
one obtains a nonlinear shallow water approximation of the energy equation. Arc 
these two approaches equivalent? We show that the answer is no. 

First we present the energy equation in three different forms: full water wave 
equations, dispersive SWE and non-dispersive SWE. Surprisingly, the energy equa- 
tion is the same for dispersive and non-dispersive SWE at leading order. Then we 
present some numerical computations over a flat bottom. It allows us to concentrate 
on the generation process and the energy transfer through the moving seabed. We 
refer to Duytkh ct al. [2008] for simulations of some real world events including 
energy pumping and its transformation over uneven bathymetry, using the SWE 
with energy. Finally we present the energy equation in a fourth form: the linearized 
dispersive water wave equations. We solve the Cauchy-Poisson problem of tsunami 
generation. Exchanges between potential and kinetic energies are clearly revealed. 



2. Derivation of the energy equation 

Consider the 3D fluid domain shown in Figure 1. It is bounded above by the free 
surface z* = rj*(x* , y*,t*) and below by the solid boundary with prescribed motion 
z* = —h*(x*,y*,t*). A Cartesian coordinate system with the z*— axis pointing 
vertically upwards and the x*Oy*— plane coinciding with the still- water level is 
chosen. 

The fluid is assumed to be inviscid. Its motion is governed by the 3D Eulcr 
equations, written here in their incompressible form (see for example Gisler [2008] 
for the compressible counterpart, after replacing p by p in the pressure term of his 
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Figure 1. Sketch of the 3D fluid domain for wave generation by a moving bottom. 
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where p* is the fluid density, u* = (u* ,v*,w*) the velocity vector, e* the sum of 
kinetic energy density e* K = | |u*| 2 and potential energy density e* p = gz*, p* the 
pressure and g the acceleration due to gravity. In the present study g = (0, 0, —g). 
For incompressible flows, the energy equation (2.3) is redundant. Indeed it can be 
obtained from Eq. (2.2). However we keep it since it is not equivalent to derive 
shallow water equations with or without the energy equation. 

Equations (2.1)-(2.3) have to be completed by the kinematic and dynamic 
boundary conditions. Since surface tension effects are not important for long waves, 
the dynamic boundary condition on the free surface reads 



P 



Psi a t z * 



V 



(2.4) 



Later we will replace the surface pressure p* by but we keep it arbitrary for now. 

The kinematic boundary conditions on the free surface and at the seabed are, 
respectively, 

* dif „ di]* drj* „ „ 

w = — h u — h v — — , z = ri , (2.5 

dt* dx* dy* ' y ' 

dh* ,dh* ,dh* 

w ' = —-^ u v — ) 

at* ox* ay* 

Below we denote the horizontal gradient by and the horizontal velocity by u^. 
After a few manipulations and integration across the water column from bottom to 



= -h* 



(2.6) 
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top, one can write the following global energy equation: 

dE 



dt* 



+ V_l • $* + P* = 0, (2.7) 



where E* is the sum of kinetic and potential energies in the flow, per unit horizontal 
area, the horizontal energy flux vector, and P* the net energy input due to the 
pressure forces doing work on the upper and lower boundaries of the fluid. They 
are given by the following expressions: 

E* = £ (lp*\u*\ 2 + p*gz^dz*, (2.8) 

$* = y" ul (^p*\u*\ 2 +p*+p*gz^dz*, (2.9) 

P* = P* s r,* t ,+plK„ (2.10) 

where p* is the pressure exerted on the free surface and pi the bottom pressure. In 
the case of a stationary bottom boundary and of a free surface on which the pressure 
vanishes, then as expected the net energy input P* is identically zero. Energy can 
be brought to the system by a moving bottom or by a pressure disturbance on the 
free surface. From now on, we take p* = 0. 

(a) Dimensionless equations 

The problem of tsunami propagation possesses two characteristic length scales: 
the average water depth ho for the vertical dimension and a typical wavelength I for 
the horizontal dimensions. It is classical to introduce the following dimensionless 
variables. The scaling for the independent variables is 

x* y* z* Vgho,* 

In order to introduce the dimensionless dependent variables we need one more 
parameter, the typical wave amplitude a: 

h ho ho h V* , h* 

w = — w , f] = — , h — — , 



a^/gho aygho I «V gho a h ' 

p* + p*gz* p* e* 

7r = , p= — , e = — — . 

p*ga p*gh gh 

The hydrostatic pressure —p*gz* has been incorporated into ir. 

The following dimensionless parameters, which are assumed to be small, are 
introduced: 

e :— a/ho, p := ho/l. 

The parameter e represents the relative importance of nonlinear terms and p, mea- 
sures dispersive effects. Note that 

1 1 F 2 

2/ 2 i 2\ , L fc 2 i 

e = -e (tr + v z ) + --^-ur + z. 

2 2 p z 
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The Euler equations of motion (2.1)-(2.3) become in dimcnsionless form 

/J?(u x + v y )+w z = 0, (2.11) 

\j?u t + en 2 (u 2 ) x + en 2 {uv) y + e(uw) z + p?ix x = 0, (2-12) 

l?v t + e/j?(uv) x + ep 2 (v 2 ) y + e(vw) z + i?-k v = 0, (2-13) 

\j?w t + sn 2 (uw) x + ey 2 (vw) y + e(w 2 ) z + [i 2 ir z = 0, (2-14) 

H 2 e t +en 2 ((e+p)u) x +en 2 ((e + p)v) y + e ((e + p)w) z = 0. (2.15) 
In dimensionless form the boundary conditions (2.4)-(2.6) become 

7T = 1], z = er], (2-16) 

w = [i 2 r\ t +efi 2 ur] x + e{i 2 vi] y , z = er), (2-17) 

ew = —]i 2 h t — e[i 2 uh x — e]i 2 vh y , z = —h. (2.18) 

In the case of a static bottom h = h(x, y), the order of magnitude of the vertical 
velocity w at the bottom is 0(/j, 2 ). With a moving bathymetry the behaviour is 
different: 

w\ z= _ h = -^h t + 0(» 2 ). 

(b) Integration over the depth 

Next we reduce the above 3D problem (2.11)— (2.18) into a 2D one by integrating 
all the equations over the water column. We begin with the continuity equation 
(2.11) which we integrate with respect to z from —h to en. Taking into account the 
boundary conditions (2.17) and (2.18), one obtains 

erj erj 

Vt + faj udz+ dy J V dZ = ~\ ht ' ^ 2 ' 19 ^ 

-h ' -h 

A source term appears in equation (2.19) due to the moving bathymetry. For 
tsunami generation, the function h* can be represented as follows: 

h*(x*,y*,t*) = h*(x*,y*) - C{x*,y*,t*), 

where h^(x* , y*) is the static sea bed profile and (*(x* , y* , t*) the bottom displace- 
ment due for example to coseismic displacements or landslides. In nondimensional 
form this representation takes the form: 

h(x, y, t) = h (x, y) - e((x, y, t), 

since the bottom variation must be of the same order of magnitude as the typical 
wave amplitude a: (*(x* ,y* ,t*~) = a((x,y,t). Differentiating with respect to time 
yields 



e dt dt y ' 



Below we will replace -h t by — ( t 
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Integrating the vertical momentum conservation equation (2.14) yields 

erj erj er\ 

d f d f d f 

*\z=-h = 7 l+Q- t J WdZ + £ fa J UW dZ + £ dy~ J VWdZ ' (2 ' 20) 

— h —h —h 

A more general expression for the pressure can be obtained if we integrate 
equation (2.14) from z to er) and use the boundary conditions (2.16) and (2.17): 

erj erj erj 

d f d f , d 



7r = t? + — / w dz + e— I uw dz + e— / vw dz -w . (2.21) 

dt J dx J dy J n 2 

z z z 

The vertical velocity w is obtained by integrating the continuity equation (2.11) 
from —h to z and applying the seabed kinematic condition (2.18): 

w = At - ^ J udz + |- j 'v dz) . (2.22) 

V -h -h / 

Finally, the integration of the equation of energy conservation (2.15) yields 

sr/ erj erj 

dtj edz + £ faJ ( e +P) udz + £ Q^ J (e+p)vdz-ep\ z= _ h ( t = 0, (2.23) 

— h —h —h 

which is nothing else than the dimensionless counterpart of Eq. (2.7). 

All the equations derived above are exact and no assumption has been made 
about the orders of magnitude of e and fi. 

(c) The nonlinear shallow-water equations with energy equation 

In Appendix A, we briefly summarize the derivation of various systems of 
shallow-water wave equations. The non-dispersive SWE are obtained by taking 
the limit fi — > 0. Here we provide the dispersive and non-dispersive SWE in their 
conservative forms [Eskilsson and Sherwin, 2005] with dimensions, based on the 
depth-averaged horizontal velocity u* . The total water depth h* + n* is denoted by 
H*(x*,y*,t*). The definition of E* is given by Eq. (2.8). 

(i) SWE with dispersion and energy equation 



dH 
~dt 



+ V • (W ®u* + 1 -gH* 2 j 



6 V V h* 
dE* 



,*2 



dt* 



V • u* IE 



1 



p*gH* 



(2.24) 



V(V-(fl"*u*)) = gH*Vh*, (2.25) 



P*9H*%. (2.26) 
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(ii) S WE without dispersion and with energy equation 



^-+V-(iru*) = 0, (2.27) 
d{H*u*) +v /^* Q * 0Q * + l gH * 2 ] = gH*Vh*, (2.28) 



dt* V T 

^ + V.(u'(f + ^« 2 )) = P*gH*^. (2.29) 

Between (i) and (ii) , only the equation for the evolution of the horizontal velocity 

(2.28) differs. In particular, it is interesting to note (this is the first main result 
of the present paper) that the energy equations are the same in the dispersive 
and non-dispersive cases. Differences appear only at higher order (terms of order 
0(e/i 2 )). The physical meaning is that the vertical velocity as well as the non- 
uniform structure of the horizontal velocity appear only at the next order. An 
additional remark can be made about the hyperbolic structure of the system (2.27)- 

(2.29) . We restrict our observations to one space dimension. Let us assume that 
a shock wave propagates from left to right at velocity s > 0. The states before 
and after the discontinuity are denoted by (H* r , u ; * r , E* r ) correspondingly. From 
entropy considerations one can conclude that for an admissible state > H*. 
For a general system of conservation laws v t + d x f(v) = 0, the Rankinc-Hugoniot 
relations have the simple form f(v r ) — f(vi) = s(v r — vi). 

Some simple algebraic calculations yield the following relations: 



H* H* 2 P y \H* H*J u*-u* K ' 

These formulas relate left and right states connected in the (H*,u*, E*) space by a 
shock wave. The first two relations are well known. What is new to our knowledge is 
the formula (2.30) which gives an insight into the energy states in a shock wave. In 
practice, they can be used for the theoretical analysis of bores and as a validation 
test for nonlinear SWE codes (with energy equation). 



3. Simulations of energy 

Next we illustrate the main features of energy evolution in tsunami generation. 
The importance of dispersive effects strongly depends on the extent of the source 
area (the smaller the source the stronger the dispersive effects) and the ocean 
depth in the source area [Kervella et al., 2007]. We restrict our study to the non- 
dispersive SWE (2.27)-(2.29). We solve these equations numerically with a finite 
volume method [Duytkh et al., 2008]. 

While the common practice in modeling tsunami generation consists in trans- 
lating the initial sea bottom deformation to the water surface, thus neglecting all 
dynamical effects, we prefer to include some dynamics in the process in an effort to 
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Parameter 


Value 


Dip angle, S 


13° 


Slip angle, 9 


90° 


Fault length, L* 


18 km 


Fault width, W* 


14 km 


Fault depth 


5 km 


Slip along the fault 


10 m 


Poisson ratio 


0.27 


Young modulus 


9.5 x 10 9 Pa 


Acceleration due to gravity, g 


9.81 m/s 2 


Water depth, ho 


1 km 


Characteristic rise time, to 


8 s 


a* = log(3)/t5 


0.1373 s- 1 



Table 1. Values of physical parameters used for the energy density computations. 



be closer to what happens in reality [Dutykh et al. , 2006] . We construct the bottom 
motion by multiplying Okada's static solution (ok( x *> 2/*)t by a function of time 
[Hammack, 1973]: 

h*(x*,y*,t*) = h* (x*,y*) (l-e- a ^)Co K (^,yl- 

The parameter a* is related to the characteristic time ig under consideration. We 
chose 

1 - e a f ° = - & a* = —5-. 

3 t* 

The various parameters used in the computations are given in Table 1. 

Other time laws are possible and we refer to Dutykh and Dias [2007] for more 
details. In the present numerical computations we chose /ig(a;*,y*) = ho = const. 
This choice is not only made for the sake of simplicity. Another reason is that 
Okada's solution is derived within the assumption of an elastic half-space which 
does not take into account the bathymetry. In order to be coherent with this solution 
we assume the bottom to be flat before deformation. 

With our definition of potential energy, the total energy is nonzero both at time 
t* = and at time t* — > oo. For practical applications it is important to isolate the 
energy available to the tsunami wave. One possibility is to define the wave energy 
as follows: 

Kave = E* + \ P *gh* 2 . (3.1) 

Clearly E^ avc — both at time t* = and at time t* — » oo when the wave has left 
the generation region. From Eq. (2.29), one finds that the energy equation satisfied 

b y Kavc is 

+ V • (u* (iC ve + \p*9H* 2 \p*9h* 2 ^ = P*g(H* h*)^- (3.2) 

t Okada's solution is a steady analytical solution for the seafloor displacement following an 
underwater earthquake, based on dislocation theory in an clastic half-space [Okada, 1985]. 
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(a) Free surface (b) Total energy (c) Potential energy 



Figure 2. Tsunami generation leading to a dipolar wave form; t* — 4 s. The same scales 
are used in Figures 2-9: for the free surface, red denotes a rise of sea level; for the total 
energy, the scale goes from (blue) to 5 x 10 4 (red); for the potential energy, the scale 
goes from (blue) to 10 4 (red). 





(a) Free surface (b) Total energy (c) Potential energy 

Figure 4. t* = 10 s 

Figures 2-9 show the distributions of free-surface elevation r/*(x* 7 y*), total wave en- 
ergy E^ ave (x*, y*) and potential wave energy \p*grf 2 at various times. One clearly 
sees the generation process. The formation of the leading elevation and depres- 
sion waves takes a few seconds. Then the propagation begins. As shown by Ben- 
Menahem and Rosenman [1972], tsunami energy radiates primarily at right angles 
to a rupturing fault (sec also Kajiura [1970]). The distribution of potential energy 
makes sense when one compares the potential energy plots with the free-surface 
plots. The total energy spreads in a more uniform manner across the area affected 
by the waves. 
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(a) Free surface 



(b) Total energy 
Figure 5. t* = 20 s 



(c) Potential energy 




(a) Free surface 



(b) Total energy 
Figure 6. t* = 40 s 



(c) Potential energy 




(a) Free surface 



(b) Total energy 
Figure 7. t* = 80 s 



(c) Potential energy 




(a) Free surface 



(b) Total energy 
Figure 8. t* = 120 s 



(c) Potential energy 
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(a) Free surface (b) Total energy (c) Potential energy 

Figure 9. t* = 160 s 




t* s t", s 

(a) Total energy (b) Total, potential and kinetic energies 



Figure 10. (a) Total energy Jj £^ avc dx* as a function of time computed with the SWE 
with energy (2.27)-(2.29) (top curve), the linearized water wave equations (middle curve) 
and reconstructed from the flow variables ti* and rj* when no energy equation is used 
(bottom curve), (b) Partition between potential and kinetic energies as a function of time 
for the reconstructed energy. The top curve of (b) is the same as the bottom curve of (a). 



It is also of interest to see how the wave energy £* ave , integrated over the 
whole fluid domain, varies during the generation process. This evolution is shown 
in figure 10. Three curves are plotted in figure 10(a). The solid curve simply is the 
energy (2.8) integrated over the water surface. The middle curve will be explained 
in Section 4 (linearized theory). The bottom curve is obtained as follows. Imagine 
that one is solving the SWE without the energy equation. In order to look at the 
energy, a natural way to do it is to approximate the wave energy E^ avc as 

Z7 1 * ^ J * r_r* I,-;* |2 i * 
^wavc ~ 2 P ' ' + 2 P9V ' 

and then to integrate over the whole fluid domain. The second main result of the 
present paper is that the difference between the proper way to compute the energy 
and the approximate way can be large. It is probably due to the vertical velocity, 
which is completely neglected in the second approach. Once the motion of the sea 
bottom has stopped, the total energy remains constant. However the reconstructed 
energy decreases with time. This indicates that computing the wave energy directly 
from the conservative energy equation (3.2) is much better, especially with the type 
of numerical method we use. In figure 10(b), one clearly sees the exchange between 
potential and kinetic energies until equipartition is reached. 
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4. Energy in the framework of the dispersive linearized 

equations 

In the case of tsunamis generated by earthquakes, nonlinear effects are not im- 
portant during the process of generation and propagation. This is why it is valid 
to use the linearized water-wave equations. Dutykh et al. [2006] and others showed 
that taking an instantaneous seabed deformation is not equivalent to instanta- 
neously transferring the seabed deformation to the ocean surface, except in the 
framework of the linearized shallow water equations (very long waves). The differ- 
ence comes from the vertical velocities and dispersion. In this case we must go back 
to the initial set of equations (2.1)-(2.3). Since the motion starts from the state of 
rest, it can be considered as irrotational (potential flow) and one can introduce the 
velocity potential u* = V0*. 

We perform the linearization of the equations (2.1) and (2.2), and of the bound- 
ary conditions (2.4)-(2.6). It is equivalent to taking the limit of the equations as 
e — > 0. For the sake of convenience, we switch back to the physical variables. The 
linearized problem in dimensional variables reads [Dutykh and Dias, 2007] 

A ± </>* +d 2 (t)*/dz* 2 = 0, (x*,y*,z*) e K 2 x [-/)*, 0], (4.1) 
dn* 

-q^i z * = 0) kinematic condition, (4.2) 

0, z* = 0, dynamic condition. (4.3) 

Within linear theory the forces that cause perturbations are so weak that the 
boundary condition at the bottom (2.6) is also simplified: 

dz* ~ dt* ' 

The bottom motion appears in the right-hand side of Eq. (4.4). 

The Laplace equation (4.1) together with the boundary conditions (4.2), (4.3) 
and (4.4) determine the boundary-value problem for the velocity potential qb* within 
the linear theory. In order to solve the equations for a prescribed bottom motion, one 
can use Fourier and Laplace transforms (this approach is followed here) or Green's 
functions. Three scenarios are considered for the bottom motion [Dutykh et al., 
2006] : the passive generation in which the deformation of the sea bottom is simply 
translated to the free surface (one is then solving an initial value problem) and two 
dynamical processes (*( x * ,U* ,t*) = T(t*)Q) K (x* , y*), where Cq K {x* ,y*) is given 
by Okada's solution. The two choices for T are the instantaneous deformation with 
T(t*) = Ti(t*) = H(t*), where H(t*) denotes the Heaviside step function, and the 
exponential law used above in Section 3: 

n*- )= ={;•_,-«,... f ;;>o! ^«*>°- 

Let C > q K be the Fourier transform of Q) K , ^ = I tanh(|k*|/i*) the dispersion 
relation and x* = (x*,y*). We provide the general integral solution for the free 
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Parameter 


Value 


Dip angle, S 


13° 


Slip angle, 9 


90° 


Fault length, L* 


150 km 


Fault width, W* 


50 km 


Fault depth 


35 km 


Slip along the fault 


15 m 


Poisson ratio 


0.27 


Young modulus 


9.5 x 10 y Pa 


Acceleration due to gravity, g 


9.81 m/s 2 


Water depth, h* 


4 km 


Characteristic rise time, £*, 


50 s 


a* = log(3)A5 


0.0220 s" 1 



Table 2. Values of physical parameters used for the Cauchy-Poisson analysis of tsunami 

generation. 



surface elevation in the three cases, 

1 



<(x*,i*) = 



Cq K e lk * x * cos(wi*) dk*, (passive) 



(4.5) 



R 2 



(2 



i_ ffChK^L 

tt) 2 J J cosh(|k*|/i* 



»?:(x*,t*) 



(2 



c, 



) 



cos(wt*) dk*, (instantaneous) (4.6) 



cos(wt*) — w/a* sin(w£*) 



cosh(|k*|/i*) 



a* 2 + cj 2 



dk*, (4.7) 



and the velocity potential in the passive case (the expressions for the other two 
active cases are a bit cumbersome) [Dutykh et al., 2006], 



S 2 



Then one can easily compute the kinetic and potential energies: 



E* K - - 2 P* 



|V0*| 2 d X *dz* 7 E* P = \p*g 
■h* 1 



rf 1 rfx*. 



(4.8) 



(4.9) 



Results are shown in figures 11 and 12. Even though there are differences during 
the first few seconds, the three mecanisms lead to the same almost exact equiparti- 
tion between kinetic and potential energy once the dipolar waves start to propagate. 
The simplest estimate proposed for the energy of tsunamis generated by a disloca- 
tion source is that given by Okal and Synolakis [2003] . They compute the increase 
in potential energy of the ocean by displacing a volume of water S x Sh* from the 
bottom to the surface of the ocean. This also represents the work of the pressure 
forces displacing the ocean bottom. Then they explain that the center of mass of 
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°0 1 00 200 300 400 500 600 700 800 °0 1 00 200 300 400 < 500 600 700 800 

Time t , s Time t , s 



(a) Kinetic energy (b) Potential energy 

Figure 11. Time evolution of kinetic and potential energies (4.9) for three mechanisms 
of tsunami generation: passive generation, instantaneous bottom motion and exponential 
bottom motion. The square in plot (b) indicates the estimate (4.11). 
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■■■ Instantaneous 
-Exponential, t = 50 s 
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(a) Total energy 



Kinetic energy E K , J ) 

(b) Trajectory in energy space 



Figure 12. Same as figure 11. Time evolution of total energy -Ewavc = E* K + E* P (a) and 
trajectory in energy space (b). The straight dashed-dotted line represents equipartition of 
kinetic and potential energies. The squares represent the last computed points. 

the displaced water, initially at height Sh* /2 above the ocean floor, is transferred 
to the ocean surface, so that the change in potential energy is not as much. The 
difference between the two is the energy available to the tsunami wave: 



// 



E* 

wave 



dx* = -p*gS(5h*y. 



(4.10) 



Incidentally, this expression does not depend on the sign of Sh* and is also valid for 
a sudden subsidence of a section of the ocean floor. It can be extended to a more 
realistic sea floor deformation, such as the one used in this paper, Q) K {x* , V*)'- 



e: 



(4.11) 



This quantity corresponds to the square in figure 11(b). 

5. Concluding remarks 

In this article we provided a formal derivation of the total energy equation in 
the framework of the nonlinear SWE, both for dispersive and non-dispersive waves. 
We also made an attempt to better understand the energy transfer from a moving 
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bottom to the water above. The importance of this topic is clear, given the serious 
hazard that tsunamis represent for coastal regions. 

Tsunami energy can be studied at several levels. A simple formula was given by 
Okal and Synolakis [2003]. In the present paper, we extended it to more spatially 
realistic sea floor deformations. But this formula does not involve any dynamics, 
which can play an important role in tsunami generation. A somewhat counterin- 
tuitive consequence of this simple estimate is that the energy carried by the full 
tsunami wave is practically independent of depth [Okal and Synolakis, 2004]. One 
would think that the energy involved in lifting a 4 km column of water above a 
rupture is different from lifting only 100 m of water (whether instantaneously or in 
a few seconds). But no difference can be found in wave height. The assumption of 
incompressibility may in fact no longer be valid, especially for a very deep ocean. 

The emphasis of this paper has been on tsunami generation. At the runup or in- 
undation stage, energy is also quite important. Potential energy is being transferred 
into kinetic energy and the study of these exchanges is left for future work. Our 
approach can also be used to analyze the structure of the wave field in caustics more 
accurately than with ray theory [Berry, 2007] . The present paper can be considered 
as a first step towards a better understanding of tsunami wave energy, in order to 
provide scales for tsunami magnitudes for example. More profound mathematical 
and physical analysis is needed. 

Another extension is the study of other mechanisms. For example Okal [2003] 
showed that the combination of a lesser absolute level of excitation and a more 
pronounced shift of the spectral energy towards higher frequencies characterized by 
strong dispersion makes landslide sources significantly deficient far-field tsunami 
generators, as compared to classical dislocations. 

The authors would like to thank Jean-Michel Ghidaglia who suggested the idea 
of considering Boussinesq equations with energy. We are also grateful to him for 
numerous fruitful discussions on numerical and fluid mechanics topics. The second 
author acknowledges the support from the EU project TRANSFER (Tsunami Risk 
ANd Strategies For the European Region) of the sixth Framework Programme 
under contract no. 037058. 



Appendix A. Derivation of dispersive shallow-water 
equations with variable bathymetry 

We derive the Boussinesq equations following the method used, for example, 
by Nwogu [1993], Villeneuve and Savage [1993], Yoon and Liu [1989]. A moving 
bathymetry was considered in the framework of the Boussinesq equations by Vil- 
leneuve and Savage [1993], but they dealt with a 2D problem leading to a ID system 
of equations. 

We need to know the depth dependence of the horizontal velocity ui in order 
to reduce the problem to a 2D one. We expand ui in a Taylor series in the vertical 
coordinate z about the seabed z = —h: 

(z + h) 2 

u(x,y,z,t) = u\ z= _ h + (z + h) u z \ z= _ h + u zz \ z= _ h + .... (Al) 

From now on, the horizontal velocity at the bottom is denoted by 

u fc := u±(x,y,-h,t). 
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If the flow is assumed to be irrotational, there is an additional relation which 
closes the system: 

(ui) 2 =Vw. (A 2) 

Substituting (2.22) into (A 2) clearly shows that (u±) z (x,y,—h,t) — 0(/j, 2 ). Sub- 
stituting (A 1) into (2.22) and integrating then yields 



^ 2 C t -M 2 V- {(z + h)u b ) + 0(^). 



(A3) 



The vertical velocity varies linearly with respect to z over the depth at leading 
order 0(p 2 ). 

The horizontal velocities can be found by integrating the irrotationality condi- 
tion (A 2) from —h to z: 

z 

u±=u b + J Vw dz. (A 4) 

-h 



Substituting (A 3) for w and integrating gives 



u_l = u b + n 2 (z + h) 



+ 0(^). (A 5) 



We see that the horizontal velocities vary quadratically with respect to z over the 
depth at leading order 0(fi 2 ). 

Following Ursell [1953], we introduce a number which measures the relative 
importance of nonlinear and dispersive effects in long waves: 

In order to simplify the computations we now assume that S = 0(1). Dispersive 
terms can be neglected if 5 3> 1. 

An expression for the pressure can be obtained by substituting (A 5) and (A 3) 
into (2.21), integrating and retaining leading order terms: 



V + ^-z 2 V • u bt + A(V • (hu bt ) - ( tt ) + 0(m 4 ). 



(A 6) 



The equation for the free-surface evolution is derived by substituting (A 5) into 
the depth-integrated continuity equation (2.19) and integrating: 

rjt + V • ((h + e V )u b ) - M 2 V ■ [y V (V ■ (hu b ) - ( t ) - y V (V • u 6 )] - & = 0. (A 7) 

The equation for the evolution of the horizontal velocity is obtained by substi- 
tuting (A 5) and (A 6) into (2.12) and (2.13): 



\h 2 1 
u bt + e(u b ■ V)u fc + Vrj + A* 2 [yV(V • u 6 ) - W(V • (hu b ) - C*)J = 0. 

Using the irrotationality condition this equation can be rewritten as 

u bt + |v|u b | 2 + Vr, + M 2 [yV(V • u b ) - W(V • {ha*) - Ct)] = 0. (A 8) 
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Finally we write the energy equation by substituting (A 5) and (A 6) into (2.23): 



1 i \ 

e dz + - (h + eij) 2 u b 



L Vd 



e(h + eri)Ct = 0. (A 9) 



The higher-order terms are of order 0(e]i 2 ). 

Another possibility for the choice of variables is to introduce the depth averaged 
velocity. The corresponding standard Boussinesq-type equations were obtained by 
Peregrine [1967] in the case of a fixed seabed. We extend the results to a moving 
bathymetry and add the energy equation. The main advantage of the depth aver- 
aged velocity consists in the fact that the continuity equation (or equivalently the 
equation for the free-surface elevation) is very simple and exact in this variable. 

Let us rewrite all equations in terms of the depth-averaged velocity defined by 



h + er) J 



u± dz. 



(A 10) 



The depth-integrated continuity equation (2.19) yields immediately 

Vt +V- {(h + er,)u) = 0. (All) 

In order to derive equations for the horizontal velocity and the energy we need 
a relation between u b and u. The desired relation is deduced directly from the 
definition (A 10) by substituting (A 5) in it: 

2 2 

u = u 6 - ^-W(V • (hu b ) - ( t ) + ^-/i 2 V(V • u 6 ) + 0(e 2 + ei? + M 4 )- 

Inverting the last equation yields 
-h 



u b = u + (1 



2 V(V-(/m 6 )-Ct)-yV(V-u 6 )J + 0(e 2 + e^ 2 + //). (A 12) 



Substituting the relation (A 12) into equation (A 8) gives the standard Boussi- 
nesq equations for a moving bottom: 



— -'-|2 , V7„ 2|~^ 2 — ^- 



u t + -V |u| 2 + V?7 + [T [— V(V • u) - -V(V • (hu) - (t) 



0. 



(A 13) 



The energy equation is obtained by substituting the relation (A 12) into equation 
(A 9): 



67] 

III' 



e dz + eV ■ 



e dz + -(h + erj) 2 | u 



e(h + e V )( t = 0. (A 14) 



Since the energy equation is redundant for incompressible flows, the linear dis- 
persion relation is unaffected by the inclusion of the energy equation. As is well- 
known, it can be improved by defining the horizontal velocity at an arbitrary level. 
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